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ABSTRACT 

One of the complexities in modelling integrated spectra of stellar populations is the 
effect of interacting binary stars besides type la supernovae (SNela). These include 
common envelope systems, CVs, novae, and are usually ignored in models predicting 
the chemistry and spectral absorption line strengths in galaxies. In this paper predic- 
tions of chemical yields from populations of single and binary stars are incorporated 
into a galactic chemical evolution model to explore the significance of the effects of 
these other binary yields. Effects on spectral line strengths from different progeni- 
tor channels of SNela are also explored. Small systematic effects are found when the 
yields from binaries, other than SNela, are included, for a given star formation his- 
tory. These effects are, at present, within the observational uncertainties on the line 
strengths. More serious differences can arise in considering different types of SNIa 
models, their rates and contributions. 
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1 INTRODUCTION 

Spectral absorption line strengths of galaxies can be used 
to probe their star formation histories (SFHs). These tech- 
niques are becoming increasingly refined and better un- 
derstood in terms of breaking under ly ing degenera c ies in 



age and metallicitv l|Worthey| 1 19941 ; iTrager et al 



Proctor fc Sanson] 



2002 - hereafter PS02; lTrager et al 



2000 



2005 



Kavirai et al.l 12007; Smi th et al.ll2008l ) and accurate inter- 



preta t ion of integrated light observ a tions (e.g. Proctor et al.l 



20041; iThomas. Maraston fc Kornl |2004| ; iLee et all 120071 : 



Serra fc Tragerl l2007h . Uncertainties in stellar population 



properties, by fitting observed spectral line-strengths to pre- 
dictions from single age, si ngle metallicitv stel lar population 
(SSP) models (e.g. PS02. ITrager et alj|2005h . The system 
of line strengths mo s t widely used is that of Li ck indices 
jWorthev et al.lll994l ; IWorthev fc Ottavianiiri997f l based on 
observations of stars taken at the Lick observatory. This is 
the system of line-strength definitions and calibration star 
data that is used in this paper, to link effectively to observ- 
able properties of galaxy spectra. 

In some cases more complex SFH fitting, to line 
strengths or co lours, is attempted using e.g. multi- 



models limit the effectiveness of these studies. One outstand- 
ing question that needs to be addressed is: how are the 
observed spectral line strengths in integrated stellar pop- 
ulations affected by abundance changes due to interacting 
binary stars? Most modellers of composite stellar popula- 
tions take into account the element yields from SNela ex- 
plosions from interacting bin ary stars (normally the Chan- 
drase khar mass, W7 model of lNomoto. Thielemann fc Yokoil 
1 19841 ') however, yields from other types of interacting binary 
star processes, such as various common envelope binaries 
that do not necessarily become SNela, are seldom taken into 
account. 

In most studies, unresolved stellar populations in galax- 
ies are characterised by their luminosity weighted average 
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pie SSPs (e.g . iHeavens et all |2004 iRossa et alj 120061 . 



ISerra fc Tragei 120071) or an exponentiall y decayin g SFH 



IIFritze-von Alvensleben fc Gerhardl 1 19941 ; iLi et al l |2004| ; 
iGanda et all l2007h sometimes with opt ional starbursts 
ijTemporin fc Fritze-von Alvenslebenl 120061 ). Empirical fit- 
ting of too many SSPs is undermined by the degeneracy be- 
tween parameters. Optical spectra of old, metal-poor SSPs 
can look remarkably similar to those of younger, more metal- 
rich SSPs, because both increasing age and increasing metal- 
licity depreciate the blue continuum. Additional constraints 
are provided by making the SFH self consistent with the 
feedback into the inter-stellar medium (ISM) from which 
stars are made. This is the approach of galactic chemical 
evolution (GCE) modelling, where elemental abundances are 
followed for a given SFH, in order to model composite stellar 
populations. 

Single SSP fitting is more robust than trying to fit the 
full SFH, and SSP fitting can be carried out for lower signal- 
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to- noise data (S/N > 15/A). However an SSP is a poor 
representation of the SFH of a real galaxy. On the other 
hand, GCE modelling may lead to a more realistic determi- 
nation of the SFH, but requires higher signal-to-noise data 
(S/N > 30A), depending on the complexity of the under- 
lying SFH. Therefore the two techniques provide comple- 
mentary approaches to determining when stars formed in 
galaxies. 

Full spectru m fitting seems a pr omising alternative to 
Lick index fitting l)Ocvirk et al.ll200a ) but there is currently 
no model that predicts the full (visible) spectrum of stellar 
populations with a range of non-solar ab undance ratios, even 
though progres s is being made with this l|Coelho et alJl2007t 
iLee et aDl2009h . In this present work we limit ourselves to 
fitting Lick indices. 

IZhang et al.l (2005) calculated the effects on luminosity 
weighted age and metallicity estimates, derived from spec- 
tral line strengths and colours, when binary interactions are 
included in the calculations of stellar evolution. They con- 
sidered the extreme case of 100% stars in binaries. Their 
models were based on the r apid binary star evolution code of 
iHurlev. Tout fc Pols! (|2002h . which predicts physical aspects 
of binary star components, such as their mass, size, separa- 
tion, luminosity and temperature etc. with time. They took 
into account the evolution and effects on spectral character- 
istics of the above physical effects of binary stars in a popula- 
tion, and showed what this meant for line strengths in SSPs 
of different ages and metallicities. Although they concluded 
that it is important to model binary stars, their results actu- 
ally show fairly small changes in most line strengths because 
of binary evolution, with changes in line strengths of up to a 
few times the typical state-of-the-art observational errors, as 
discussed in Section T3. II of this paper. For example, changes 
in H/3 were up to +0.3 7A for populatio ns with minus with- 
out binary interactions. iLi fc Hani (2008) find similar results. 
However, they did not model the yields from binary stars, 
nor did they consider non-solar element abundance ratios in 
their SSPs. 

|Pe Ponder fc Vanbeverenl |2004) combined initial prop- 
erties of stars with stellar evolutionary models (theirs and 
others) to generate the yields from stellar populations in- 
cluding binary types. They found that including the effects 
of binaries (apart from SNela) leads to small changes in the 
predicted yields. Other uncertainties in stellar evolution (e.g. 
mass-loss rates, massive-star remnant masses etc.) make at 
least as much difference to the chemical yields. 

In this paper we take a complementary approach and 
consider the effect of incorporating yields from single and 
binary star populations into a GCE model in order to deter- 
mine the effect of chemist r y on a bsor ption line strengths. We 
apply the code of llzzardl (|2004h and llzzard et al.1 (|2006l ) to 
calculate the yields of both single- and binary-stellar pop- 
ulations of different metallicities as a function of age. In 
that work, as expected, the main difference between single- 
and binary-star populations is the effect of SNela, but yields 
from giant stars are also reduced in binary populations. This 
is because when a star becomes a giant it is most likely to 
interact with a close companion and lose mass through non- 
conservative Roche-lobe overflow and common-envelope evo- 
lution. The yields of elements produced in asymptotic giant 
branch (AGB) stars e.g. carbon, nitrogen and s-process ele- 



ments such as barium, are reduced by up to 20% compared 
to single stars. 

We combine the single- and binary-star yields with a 
galactic chemical evolution model which parameterises star 
formation histories with six variable parameters. These are 
tuned to simulate a range of galaxy types (PS02). 

Unfortunately, a major uncertainty in our yield calcu- 
lation, and hence in our GCE models, arises from the cal- 
culation of yields from type la supernovae. While our bi- 
nary population synthesi s model predicts a certain la rate 
based on the algorithm of lHurlev. Tout fc Polsl (|2002h . there 
are several progenitor channels, all of which are seriously af- 
fected by model parameters (e.g. the common-envelope ejec- 
tion efficiency). Worse still, it is not known which progenitor 
channels are truly responsible for SNela. 

In order to consider quantitatively this uncertainty, we 
have split our binary yield set into two parts: SNela and 
"other binary yields" . The SNela set is then subdivided into 
two channels: 1) Sub-Mch explosions, which in our model 
arise from edge-lit detonations (ELDs), where Mch ~ 1.4Mq 
is the Chandrasekhar mass. These explode when a CO white 
dwarf (COWD) accretes more than O.15M0 of helium-rich 
material. 2) Mch explosions arise from the accretion of ma- 
terial onto a COWD until it reaches a mass of Mch- This 
set includes COWD-COWD mergers. 

The "other binary yields" set, which includes yields 
due to all other processes in binary stars such as non- 
conservative Roche-lobe overflow, common- envelope ejecta, 
novae, stellar winds etc., is then treated separately from the 
uncertainties in SNela yields. 

There is a time delay between the formation of a binary 
system that could lead to a SNIa and the explosion itself. 
This time delay depends on the type of SNIa p rogenitor 
(e.g. |Pe Ponder fc Vanbeverenl2004l : iGreggiollioba ) and this 
affects the timescale of enrichment of the ISM with iron in 
galaxies. Thus it is important to study the effects of different 
SNela progenitors and the parameters associated with their 
production rates. 

Previous authors have described how the various 
sources of enrichment are expected to co ntribute to the 
yields with time in a stellar population (e.g. IWorthevlll998l ; 
iBallero et al. I l2007l ). Briefly, type II SNae from massive stars 
will contribute a broad range of metals including the im- 
portant aip/ia-capture elements (0,Ne,Mg,Si,S.Ar,Ca etc.) 
on short timescales (^ 5 x 10 7 years, approximately the 
pre-supernova lifetime of an 8M0 star); whereas SNela 
from interacting binary star channels contribute mainly 
iron-peak elements (Cr,Fe,Ni), over lo nger timescales (e.g. 
M annucci, Bella Valle fc Panagialbood ). Intermediate mass 
single stars recycle H and He and also contribute to He, 
C, N and O yields. Intermediate mass stars in binary sys- 
tems other than SNIa al so contribute small amounts of 
other heavy eleme nts (e.g. |Pe Ponder fc Vanbeverenl 120041 : 
llzzard et al.ll200"?3 ). These different sources are all included 
in our GCE models. 

This paper addresses the question of how these yields 
affect spectral line strengths in self-consistent GCE mod- 
els, covering a variety of star formation histories. Variations 
of abundances and abundance ratios are taken into account 
through the SSP line-strength models used. Section [5] de- 
scribes our models of binary-star yields and our GCE mod- 
els. Section [3] describes how these combined models are used 
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to test the effects of including other binary yields, besides 
SNela, on line strengths; Section [4] explores some effects of 
different SNIa models and common envelope parameters on 
line strengths; some discussion is given in Section[5]and Sec- 
tion [6] summarises our results. 



2 MODELS 

In this section we briefly describe our models. We begin 
with our single- and binary-star stellar yield calculations 
and follow with a description of their implementation into 
our galactic chemical evolution code. 

2.1 Yields from single and binary stars 

Our yields are calculated for single and binary stars accord- 
in g to the rapid bina ry evolution and nu cleosynthesis models 
of llzzardl |2004 ) and llzzard et alj (120061). Stellar evolu tion is 
based on the algorithm of lHurlev. Tout fc Polsl l|2002l ') which 
uses analytic fits to stellar properties such as luminosity, ra- 
dius and core mass to follow stellar evolution for stars of 
masses up to IOOMq from the main sequence to the end of 
evolution, either as a white dwarf or in a supernova explo- 
sion. This is coupled to a binary interaction scheme which 
models tidal interaction, mass transfer by winds and Roche- 
lobe overflow, mass accretion, common envelope evolution, 
novae etc. 

The nucleosynthesis routine runs in parallel to the stel- 
lar evolution algorithm and follows the surface abundances 
of stars during all phases of evolution. It includes the fol- 
lowing: 

• A synthetic thermally pulsing asymptotic giant branch 
model which mimics the behaviou r of the detailed models of 
iKarakas. Lattanzio fc Polsl (|2002l ) and lKarakas fc Lattanziol 
(|2007l V including third dredge up and hot-bottom burning 
by the CNO, NeNa and MgAl cycles. 

• Surface abundances, a nd hence stel l ar wind yields, from 
massive stars according to iDrav et all l|2003al ). lDrav et al.l 
(2003b]) and Stancliffe (private communication). 

• Core-collapse supernov ae according to the calculations 
of IChieffi fc Limongil (|2004 ). 

• SNIa yields from sub-Chandrasekhar mass (sub-Mch) 
and Chandrasekhar mass (Mch) explosion s (yields from 
Livne fc Arnettlll995l and the DD2 model of llwamoto et all 
19991 respectively). 

• Nova yields of I Jose fc Hernanj (|l998T l. 

• Stellar wind collision, accretion and time-dependent 
thermohaline mixing into secondary stars. 

Many parameters are associated with both our single- 
and bi nary-star models, these are discussed in some de- 
tail in llzzard et al.l (|2006T ) . We make no changes from the 
standard assumptions used there, except that orbital en- 
ergy transfer param e ter q ce is set to 3 as suggested by 
iHuiiev. Tout fc Polsl (|2002f l in order to maximise the ef- 
fect of SNela. The common e nvelope structu r e par ame- 
ter Ace is fitted to t he data of iDewi fc Taurisl (2000) and 
iTauris fc Dew] (|200lf ) although we optionally allow for a 
fixed value of 0.5 (see Section |4, 111 . 

This constitutes our "baseline model" referred to in 
later sections. We note that the rate (and hence yield) of 



Mch SNela is considerably less than typically used in GCE 
models, but we make up for this with copious numbers of 
sub— Mch supernovae. A detailed study of the rates and 
yields of SNela is beyond the scope of this paper but a pre- 
liminary exploration is carried out in Section [4] 

Our yield sets are calculated by integrating yields for 
a pop ulation of stars using the iKroupa. Tout fc Gilmorel 
(1993) initial mass function (IMF) with masses between 0.1 
and 80Mq for single stars and binary primary stars, a bi- 
nary secondary mass distribution which is flat in mass ratio 
q such that all secondary masses (less than the primary mass 
i.e. q < 1) are equally likely, and a binary separation dis- 
tribution which is flat in log-separation (i.e. distributed as 
1/a where a is the binary separation) between 3 and 10 4 Rq. 
The yields were generated in lMyr bins for a maximum time 
baseline of 14Gyr. It proved difficult to temporally resolve 
yield output at late times when using our standard logarith- 
mic grid for the primary mass. To resolve this we applied an 
algorithm which chooses the primary masses as a function of 
the main-sequence lifetime such that yield rate at late times 
is time-resolved and smooth. 

We split our yield calculations into three sets: 

• Single stars (s), 

• Binary stars, excluding SNela (b), and 

• SNela (of all types) (SNIa). 

These contributions to yields are combined in different 
proportions in our GCE models (see Sections 12.21 and |3)| . 

2.2 Galactic chemical evolution models 

Our GCE model self-consistently follows the mass and chem- 
istry of both stars and the ISM for a given SFH (see PS02 
for a detailed description) which is described by a simple 
but flexible parameterisation. A Salpeter IMF is assumed 
for the stars, as available for the SSPs. A Schmidt law is 
assumed for the star formation rate: SFR = Cg, where g 
is the gas mass within a fixed volume and C is the effi- 
ciency of star formation. Sixteen of the most abundant ele- 
ments are followed individually in our GCE model and are 
then combined into groups when linked to SSP models. The 
chemistry and mass of stars thus generated as a function 
of time are converted into observables through the use of 
SSPs . We use line strength predictions from the SSP mod- 
els of lThomas. Maraston fc Kornl i|2004T ) - hereafter TMK04. 
These allow not only for age and metal mass fraction (Z) 
variations in stellar populations but also for dependencies 
of line strengths on element abundance ratios, characterised 
by the ratio of a-capture to iron-peak elements ([a/Fe] 
ratio). Following TMK04 we combine the abundant alpha- 
capture elements 0,Ne,Mg,Si,S,Ar,Ca plus N and Na as the 
a-elements group of elements in order to predict [a/Fe] for 
the SFHs modelled. Carbon remains separate from the a- 
elements in the models of TMK04. 

The parameters that describe the SFHs of our GCE 
models are listed in Table [T] and a schematic illustration is 

1 Where [a/Fe] = log(o/Fe) — log(«0/Fe0) , with a and Fe being 
the gas mass fractions of a-capture and Fe-peak elements from 
which the stars are made. See TMK04 and references therein for 
details of the elements included. 
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given in Fig. Al of Appendix [A] These allow for a wide 
range of possible SFHs, covering early star formation (SF), 
described by the SF efficiency (Co) and gas inflow rate (Fo), 
plus an optional, more recent SF event. There are param- 
eters for the onset time (Ti), duration (Di) and star for- 
mation efficiency (Ci) of this later event, as well as the gas 
inflow rate (Fi) during this more recent SF event. The SF 
and gas flow are switched off after the merger event be- 
cause elliptical galaxies show little evidence for ongoing star 
formation or massive gas flows. This star formation shut- 
off (often referred to as "quenching") mimics the fact that 
the remaining gas has been made unavailable for star for- 
mation throug h some process, either Active Galactic Nu- 
cleus feedback (|Bower et all l200rj ; [Hopkins et all 120071), or 
heating through a ccretion shocks l|Dekel fc Birnboimll2006l ; 
lOcvirk et alboosl ). Units for these parameters are also given 
in Tabled] 



l|Hau.Carter fc Balcellsl Il999h . The SFH models explored 
here do not attempt to model recent SF or late-type galax- 
ies. Instead we concentrate on modelling the effects of bi- 
nary star yields on galaxies that are early-types now, for 
which the question of binary star effects on absorption line 
strengths is more pressing due to the slow changes of spectral 
features with stellar population parameters. These SFHs are 
listed and described in Table [1] (and sh own later, in Fi g. [3j) . 
The binary fraction is uncertain (e.g. I Converse fc Sta hlcr 
2008 and references therein). Here, for simplicity, we as- 
sume 50% mass fraction in single stars and 50% in binary 
stars. First, the effects of yields from binary stars other than 
SNela explosions were tested (Section [3}, using normalisa- 
tions based on recent observations (described in Section l3~Tj) . 
Then the effects of different assumptions about SNela pro- 
genitors and common envelope parameters were tested (Sec- 
tion ED. 



Our GCE model incorporates tables of yields from stel- 
lar models. For the present work, the yields used are from the 
single and binary star SSP models described in Section \2. II 
The GCE code models the mass of stars and their element 
abundances for populations in each time step. Appropriate 
SSP model information (integrated line strengths and lumi- 
nosities) is chosen for the relevant age, metallicity and [a/Fe] 
ratio for stars at each time step. The luminosity weighted 
sum of SSPs then predicts the observable characteristics, i.e. 
the line strengths, of the integrated stellar population. 

There are no SSP models of line strengths that incor- 
porate both binary evolution effects and effects of non-solar 
abundan ce ratios. Th e form er effect has been shown to be 
small bv lZhang et"al] l|2005h . The latter effect is important 
for modelling accurate line strengths in galaxies. Our SSP 
line-strength models are based on TMK04 because they 
allow a variable [a/Fe]. However, TMK04 do not predict 
broadband colours, so our luminosity weightings (in B, V 
and I ) are taken from the SSP models of lBruzual fc Charlotl 
l|2003l . see|http://www.cida.ve/~bruzual/bc2003|). Their lu- 



minosities for a Salpeter IMF are used, for consistency with 
the SSP line strength models, and assuming Padova 1994 
tracks. The filters used to compute the B, V, I luminosities 
are respectively the Buser B3 filter, Buser V filter and the 
Cousins / filter. 

Spectral line strengths are thus predicted for integrated 
populations containing yields from binary and/or single 
stars. Our GCE models have a fixed time step (30Myr) 
and line strengths are calculated at 13.7Gyr - the approx- 
imate age of the U niverse from the WMAP 3 year results 
i|Spergel et al.ll2007l ). 

To test the effects on yields and line-strength from bi- 
nary stars, a range of simple SFHs were tested both with 
and without binary stars (b above). The motivation for 
our choice of SFHs was to cover a range of representa- 
tive histories approximating early-type (E,S0) gala xies (e.g. 
iRakos. Schombert fc OdenlT2008l . iReda et al.l l2007ft . includ- 
ing early (primordial) collapse models, early collapse fol- 
lowed by later star formation to mimic later mergers (at 
2 or 5 Gyrs ago) and finally a model approximating the 
rapid merger 2 Gyr ago of two spiral galaxies like the Milky 
Way in their pre-merger SFR, accompanied by a merger 
induced starburst. An example of a system which has un- 
dergone a major merger in the last 2 Gyr is NGC 2865 



3 EFFECTS OF BINARY YIELDS OTHER 
THAN SNEIA 

Canonical galactic chemical evolution models treat all stars 
as single stars and implement binary effects with ad-hoc 
prescriptions for type la supernovae and, rarely, yields from 
novacQ- Our calculation differs somewhat in that we calcu- 
late SNIa and other binary star yields (described in Sec- 
tion I2.1|l directly from our binary population simulations, 
so an ad-hoc model is not required. Therefore we implement 
true binary yields, including the SNIa yields, as calculated 
in section [2~T1 We consider a binary fraction of 50% by masfl 
such that our yield sets are: 



— s + —b + — SNIa (hereafter population B) and (1) 
1 



s + — SNIa, (hereafter population S) 



(2) 



where s, b and SNIa are our yields from single stars, 
binary stars without SNela, and SNela yields only, respec- 
tively (see section [2~Tj) . With these yield sets we simulated 
the chemical evolution of galaxies with various star forma- 
tion histories (Table [TJ to compare the effect of binary star 
yields other than SNela on fine strengths. Our GCE code 
was run with each of the above two combinations of stel- 
lar populations (given in Equations 1 and 2), for the differ- 
ent SFHs. Differences between the two assumed populations 
were thus found for gas properties and stellar line-strengths. 



3.1 Line strength changes relative to 
observational errors 

The relative differences in line strengths between popula- 
tions S and B, as a percentage of typical observational er- 
rors, are shown in Fig. [T] Table [5] shows the adopted errors 
which are described in detail in Appendix [Bl 

All five SFHs lead to similar differences in line-strength 

2 Nova yields are probably only relevant for minor isotopes such 
as 13 C. 

3 Equivalent to ~ 40% binary fraction by number of systems 
assuming our default initial distributions of binary parameters. 
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Scenario 


Co 


Fo 


Ci 


Ti 


Fi 




Comment 




Gyr" 1 


M Q Gyr" 1 


Gyr" 1 


Gyr 


M Gyr" 1 


Gyr 




PCI 


4.0 


4 x 10 6 


4.0 


0.25 





0.25 


Early, rapid SF lasting 0.5 Gyr 


PC2 


4.0 


1.3 x 10 6 


1.0 


0.75 





0.75 


Early, rapid SF lasting 1.5 Gyr 


LB1 


4.0 





4.0 


8.7 


2 x 10 6 


0.5 


Early, rapid SF + burst 5 Gyr ago 


LB2 


4.0 





4.0 


11.7 


2 x 10 6 


0.5 


Early, rapid SF + burst 2 Gyr ago 


MMW 


0.04 





4.0 


11.7 


2 x 10 6 


0.5 


Early, slow SF + burst 2 Gyr ago 



Table 1. Parameters that define the five star formation histo ries considered in this paper. The star formation 
efficiencies Co and Ci are defined as in lSansom fc Proctorl lll998l l. while Fo and Fi are infall rates in Mo /Gyr, with 
the initial gas mass being 10 6 Mq. Ti and T2(= Ti + Di) are the transition times in Gyr where the corresponding 
infall rates and star formation efficiencies apply (i.e. Ci and Fi apply during Di). The motivation for choosing these 
parameters is discussed in Appendix [X] Each set of parameters corresponds to a different evolutionary scenario. 
From top to bottom: primordial collapse with infall and star formation shutting off early (PCI) and slightly later 
(PC2). The next two scenarios describe an early- type galaxy accreting a significant amount of gas (e.g. through a 
merger) at a late (LB1) or very late (LB2) age, and where this event triggers a burst of star formation that shuts off 
0.5 Gyr after the encounter. The last scenario (MMW) represents a slowly forming, Milky Way-type disc accreting 
a large mass of gas that triggers intensive star formation at a late age. In each case the infalling gas doubles the 
mass. In these models the present age of Universe is assumed to be 13.7 Gyr. 



indices between populations B and S. The differences were 
less than 0.9 x typical observational errors, or less than 4% 
of the index strengths themselves (for indices that remain 
positive) . 

If we assume [a/Fe] = in our GCE model then the 
differences between populations B and S are smaller, less 
than 0.75 x the typical error. This reflects the fact that the 
metallicity and the [a/Fe] ratio are dominated by SNela and 
that the overall metallicity changes very little between GCE 
models which use yields B (Equation [1]) or S (Equation [2}. 

Line strengths are essentially pseudo-equivalent widths 
because a true continuum cannot be defined in complex 
galaxy spectra. For a given spectral feature, a pseudo- 
continuum i s defined by intensity levels within specified side- 
bands (e.g. IWorthev et al.l 1 19941 ) which means a given line 
strength can be positive or negative (e.g. H7 and H<5). For 
this reason the differences in line strengths presented in Fig- 
ure [T] are normalised to typical observational errors. Fig- 
ure [T] shows that the differences between line strengths cal- 
culated with populations using the yield sets B and S are 
all well within typical observational errors for all the SFHs 
tested. Fig. [2] illustrates these differences in terms of frac- 
tional changes, for those absorption line-strength that take 
only positive values. 

We note that some small, systematic effects are present 
in the line-strength differences between B and S popula- 
tions. Line strengths calculated for yield set B typically have 
slightly smaller metal-sensitive features (coloured blue, or- 
ange and black in Fig. [1} and slightly elevated hydrogen line 
strengths (coloured magenta in Fig. [TJ. Calcium (in black, 
labelled Ca) and iron sensitive features (in blue, labelled Fe) 
are depreciated by typically 0.2 x the line strength errors. 
Magnesium and carbon sensitive features (in orange and 
black respectively) show the most negative differences which 
implies that a-capture metal absorption features are weaker 
in our models that include binary yields. These changes are a 
result of enhanced mass loss in binaries which leads to ejec- 
tion of material which has been less processed by nuclear 
burning than in equivalent single stars. 

In future, and with the advent of higher signal-to-noise 
data (by at least a factor of 3), these systematic effects might 
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0.070 


TiOi (mag) 
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0.0020 



Table 2. Observati onal line-strength error s and their 

sources: PS02 = [Proctor fe Sansoml 12001 S-B0 6 = 

ISanchez-Blazquez et al.l |2006| . Dea05 = iDenicolo et al.l l2005i 
Our selected errors are in bold, as used in Fig. 1 and described 
in Appendix FBI 

become detectable, as small offsets to best fit population 
models. However, other model uncertainties currently out- 
weigh this, such as the yi e lds from massive star models (e.g . 
Wooslev fe Weaver! Il99gt |ffirschi. Mevnet fe Maederl 12003 : 
Eldridge. Izzard fe Toutl l200cf ) and the still -unknown pro- 
genitors, and associated yields, o f SNela l|lwamoto et al.l 
1 19991 : lYungelsonll2004l : iToutl I2005T I . Similarly, uncertainties 
in calibrations to the Lick standard system and in SSP 
models limit our ability to probe such small effects in real 
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populations. More accurate line-str ength systems, based on 
well ca librated stellar libraries (e.g. ISanchez-Blazguez et alj 
l2006al ). will allow line strengths to be measured more accu- 
rately and might help us to see such subtle differences in the 
future. 



3.2 Time dependent abundances 

Fig. [3] shows the differences between ISM abundances for 
populations with (B) and without (S) the binary yields other 
than SNela, for the five GCE models from Table [T] as a 
function of time. Similar levels of differences are seen in all 
the SFHs modelled with the GCE code. Abundance differ- 
ences are small (less than 7% in Z, Mg and Fe and less 
than 25 % in C and N) and follow the systematic differences 
reflected in the line strengths. That is, iron, magnesium, car- 
bon, nitrogen and the total metallicity are typically slightly 
lower in abundance in the ISM of the population including 
other binary yields. Again this is because of binary-enhanced 
mass-loss of, relative to single-stars, relatively unprocessed 
material. 

Oscillations in ISM gas abundances due to small num- 
bers of stars in the yield calculations can be seen when the 
gas mass and feedback from older stars are small. The GCE 
model assumes a gas inflow composition the same as the ex- 
isting ISM which aggravates the problem. For example gas 
inflow in the late burst (LB) models only starts at Tl, so the 
composition of the ISM then strongly influences the subse- 
quent chemistry in our models, even if there was hardly any 
ISM left just prior to Tf . This problem is drastically reduced 
(but not totally eliminated) by using sufficient numbers of 
stars in the evaluations of population yields and fine time 
steps in those evaluations when the evolution is most rapid 
(see e.g. Fig. [2]where the LB models behave overall in a sim- 
ilar way as the other models plotted there, and Fig.[3]where 
the oscillations in yields are relatively small). The problem 
does not affect the vast majority of stars in our simulated 
galaxies because the rate of star formation is proportional 
to the gas mass. 



4 DIFFERENT SNEIA 

The contributions of Fe peak elements from SNela en- 
richment are crucial to the correct understanding and in- 
terpretation of galaxy spectra through GCE modelling. 
Both the timescale and quantity of the enrichment pat- 
tern from SNela are important. Most modellers of stel- 
lar populations (integrated or resolved) ignore the uncer- 
tainties in SNela models and use the predictions of yields 
from the stand ard COWD deflagration model for SNela 
(W7 model of iNomoto. Thielemann fe Yokoil 1 1984). De- 
tailed binary-star models (e.g. IPe Ponder fe Vanbeverenl 
12004 lHan fe Podsiadlows ki 2004; Toutll2005l . and references 
therein; iMartin. Tout fe Lesa ffrc 2006) indicate how diffi- 
cult it is to form SNela through this conventionally assumed 
route of a COWD increasing its mass through steady ac- 
cretion and burning of material from an interacting com- 
panion star until it reaches the Chandrasekhar mass limit 
(Mch — 1.4 Mq). The difficulty arises in the unstable na- 
ture of the accretion: too rapid and a common envelope may 
form, which does not necessarily lead to a SNIa; too slow and 



novae (probably) lead to mass loss. The theoretical mass ac- 
cretion rates under which a standard SNIa can form cover 
only a factor of a few up to a few x 10 -7 Mq yr -1 (|ToutJl2005l 
and references therein) . 

The results of Hachisu. Kato fe Nomotol |l996h and 
lHachisu. Kato fe Nomotol (|2008) suggest that accretion 
rates of ~ 10" 7 to 1O" 5 M yr" 1 could lead to SNela if the 
accreting white dwarf has a strong wind that effectively re- 
duces the accretion rate to that required for steady accre- 
tion. Here we refer to this as the 'disk wind' model. 

There are other channels that may lead to a 
Chandrasekhar-mass COWP and hence a la supernova. 
These include close double-degenerate (COWP-COWD) 
systems, which merge because of orbital shrinkage caused 
by gravitational radiation. Theoretically, He white dwarfs in 
binaries may also contribute to SNIa, however, evidence fo r 
their existence as SNela progenitors is lacking (|Toug[2005). 

Sub- Chandrasekhar mass models, while currently out 
of favour as la progenitors, dominate the rate of SNela in 
our binary models. These are mostly edge-lit detonations 
(ELPs) of COWPs accreting 0.15M Q of helium-rich mate- 
rial which ignites, setting off the whole star. If these types of 
SNela really do dominate the yields then it is important to 
study whether there are any chemical signatures that could 
be used to confirm their presence in integrated stellar pop- 
ulations. 

Table [3] shows the relative contributions from different 
types of SNela from our binary yields at different metallici- 
ties assuming our baseline binary-star plus single-star model 
("Section 12. 1[) . This shows that yields from ELPs dominate 
the SNIa contributions in our models, particularly at high 
metallicity. The exact mix of SNIa types does not matter for 
the tests of the impact of other binary yields (from common- 
envelope evolution, novae etc.) on line strengths (Section l3.ll 
above), but it may limit our ability to accurately interpret 
line strengths from galaxies. Therefore we go on to try to 
test the impact of different SNela on spectral line strengths. 

As expected, iron dominates the yields from the SNela 
in these models (see Table [3} . Fig. 2] shows the time evolu- 
tion of iron for two different categories of SNela. Note from 
these plots how the distributions are sensitive to the initial 
metallicity of the stars. For all the possible SNela model 
yields shown in Fig. [4] the feedback from SNela includes a 
large peak at very early times (~ 10 s yr) plus extended feed- 
back at a lower level, generally decreasing over the rest of 
the time. This is in contrast to what was often assumed until 
recently, where the SNela contributions were not expected 
to contribute significantly in galaxies until a few x 10 8 yr 
IIMatteucci fe Recchj|200lh . It is, however, much more in 
agreement with rece n t ideas about SNela time s cales (e.g. 
lAubourg et~al1 |2008l; lHachisu. Kato fe Nomotol 120081 and 
references therein; Kobavashi fe Nomotol 120081 ). in which 
SNIa occur with short and long delay times. 



4.1 The effect of different SNIa progenitors and 
common envelope parameters 

In this section we carry out a series of tests to examine 
the effects on line strengths from different combinations 
of SNIa types and common envelope parameters using 
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Figure 1. Differences in indices found for populations with- without other binary yields (apart from SNela explo- 
sions), for the five different SFHs, plotted in the order given in Table [T] The differences are plotted as a percentage 
of typical observational errors. The indices plotted include the line strengths and molecular band features listed 
in Table [2] Hydrogen features are labelled in magenta (starred points), iron sensitive features are in blue (filled 
diamond points), magnesium sensitive features are in orange (pluses) and other features are in black (pluses). 



each the SFHs given in Table [T] Our GCE model was 
modified to allow yields from different proportions of 
the various SNIa types. Yield sets were calculated from 
weighted sums of the following subsets of the baseline yields: 

(i) Single stars, 



(ii) Binary stars (excluding SNela), 

(iii) SNela (baseline model), 

(iv) He WD SNela (accreting and merging), 

(v) ELD SNela, 

(vi) Others (COWDs which accrete or merge such that 
Mwd > M Ch ) 
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Figure 2. Same as for Fig. [T] except that here the fractional differences are plotted for those line-strengths that 
take only positive values (i.e. excluding H7, H<5 and CN bands). 



where (iii)=(iv)+(v) + (vi). 

That is, the sum of equal weight contributions from iv, 
v and vi above add up to the baseline model SNIa contri- 
butions (iii above). Option vi labelled 'Others' consists of 
yields from Chandrasekhar-mass SNIa explosions. 

Table [4] lists the eight test cases with corresponding 
weighting factors. The results of each test are then com- 
pared with the baseline model (which uses proportions given 
in Equation [1] equivalent to test 1). Differences that arise 



are then due to the proportions of SNIa types included. The 
maximum difference in any line strength, relative to the rel- 
evant observational error and the root mean square differ- 
ence in all line strengths relative to observation errors are 
also shown in Table [4] 

Tests 1 and 2 are equivalent to our populations B, the 
"baseline model" , and S (Eqs. Q] and [2] respectively) . Again 
this shows that the effect of changing from a purely single- 
star population with SNela to a true half-single and half- 
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Figure 3. Normalised SFH, normalised gas mass, total metallicity in the ISM and differences in ISM abundances 
for various elements as a function of time in Gyr, for the five star formation history models listed in Table [T] Upper 
row shows results for SFH models PCI (left) and PC2 (right), central row shows LB1 (left) and LB2 (right) models 
and lower row shows MMW model. The SFH axis is the rate of mass going into stars per Gyr, normalised by the 
initial mass of gas. The gas mass fraction (F 9 ) is also normalised by the initial mass of gas. Total metallicities are 
in mass fractions and abundance differences are in percentage, i.e. 100(B-S)/B, where B and S are the populations 
with and without other binaries respectively (see Equations 1 and 2). N.B. The oscillations in our models with late 
bursts (LB1 and LB2) are a result of low gas masses at those times. 
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0.0029 
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0.0038 
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0.0019 


0.0016 


0.0015 


0.0014 


ELD (%) 
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29.1 
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63.3 


65.8 


Hela (%) 


66.6 


64.0 


63.2 


37.1 


28.3 
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21.8 


Others (%) 


3.3 


6.8 


6.8 


9.0 


10.9 


11.8 


12.4 



Table 3. Mass contributions of various types of SNIa progenitors to the total SNIa yields as a function of metallicity 
Z for our baseline model. The "Total" is the mass fraction returned from SNela relative to the total initial mass of 
the stellar population. Similarly, the mass fraction returned from iron is shown as "Fe" . The final three rows show 
percentages of ejecta from edge-lit detonations (sub-Mch); He white dwarf explosions and other (Mch) explosions 
e.g. accreting COWD and COWD-COWD mergers. 
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Figure 4. Time evolution of iron yields from two different groups of SNela, from binary star model yields. The 
models shown are for ELD and Other SNela as listed in Table [3] 



binary population is within the observational uncertainties 
on the line strengths. 



Parameters that quantify the behaviour of a common 
envelope (CE) formed around a close binary can affect the 
number and type of SNela that occur, particularly the num- 
ber of double degenerates that can merge to form SNela. 
The efficiency of orbital energy transfer to the CE is pa- 
rameterised by acE (see Hurley et al. 2002). The CE may 
be completely removed by orbital energy transfer to the en- 
velope. Observational evidence for this is in the existence 
of close double degenerate binaries such as 4U1820-30 or 
B2303+46 (Church et al. 2006) which are thought to have 
achieved their close orbits via previous orbital energy loss to 
a CE. Estimates of «ce vary from less than 1 up to about 3, 
with acE > 1 possible if energy sources other than orbital 
energy are involved (Hurley et al. 2002). If cvce is large then 
energy is more easily transferred to the CE until it becomes 
unbound. Conversely, if the transfer of orbital energy to the 
CE is less efficient, then the CE may be bound for longer, 
giving the stellar orbits time to decay by various means, 
which is more likely to lead to coalescence (Church et al. 
2006). However, if orbital energy transfer is too inefficient 
this essentially removes one of the potential orbital decay 



mechanisms for the binary, leading to fewer SNela via coa- 
lescence. 

Another parameter relating to CE evolution is Ace, 
a dimensionless factor in versely proportional to the initial 
binding energy of the CE (|Hurlev, Tout fe Polsl l|2002l )) . This 
parameter is important since it affects the amount of energy 
that can be transferred to the CE before it becomes un- 
bound. This affects the closeness of the two stellar cores 
remaining after the CE is lost. A value of Ace = 0.5 was 
used in Hurley et al. and we test this value compared to our 
baseline model described in Section \2. II 

Tests 3 and 4 show that the effect on line strengths of 
changing the common envelope parameters by setting either 
qce = 1 or Ace = 0.5 is about one third of the observa- 
tional uncertainty and at most of the same order as the ob- 
servational uncertainty. Therefore uncertainties in the values 
of common-envelope parameters, and their effects on SNela 
yields, can normally be neglected in GCE models of line- 
strengths in integrated stellar populations. 

We now ask if there is a significant difference between 
using the baseline model yields (dominated by sub-Mch 
ELDs) and yields from Mch SNela scaled up to eject the 
same total mass (at Z = 0.02). The results are shown in 
test 5. The Mch yields differ in both composition and rate 
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Test Scenario Proportions Max diff RMS Comment 

Number s b All SNIa ELD Others /errors diff 













only 


only 








1 


Baseline 


0.5 


0.5 


0.5 














(B, Eq. 1) 


2 


Single stars + SNIa 


1.0 





0.5 








0.9 


0.50 


(S in Eq. 2) 


3 


«CE = 1-0 


0.5 


0.5 


0.5 








1.0 


0.39 




1 


Ace = 0.5 


0.5 


0.5 


0.5 








0.65 


0.30 




5 


Standard SNIa 


0.5 


0.5 








4.237 


1.0 


0.46 


M Ch COWDs 


6 


Exclude all SNIa 


0.5 


0.5 











3.2 


1.71 




7 


2x baseline SNIa 


0.5 


0.5 


1.0 








2.0 


0.95 




8 


Disk wind 


0.5 


0.5 


0.5 








4.0 


1.97 





Table 4. Proportions used in tests and resultant differences in line-strengths compared to the baseline model, which 
includes other binary yields as well as those from SNela (see Section 12.10 . Maximum and RMS differences are given, 
expressed as a fraction of the typical observational errors on line-strengths. From Section |2. II s are single stars, b 
are binaries excluding SNela. Tests 3, 4 and 8 used different SNIa calculations. All differences are compared with the 
baseline model. See section [4. II for discussion of the sense of the differences for particular groups of line-strengths. 



from the sub-Mch yields, yet the maximum difference be- 
tween the line strengths in both models was only of the or- 
der of the observational uncertainties with RMS differences 
of only half the observational uncertainty. Therefore it is dif- 
ficult to distinguish between Mch and sub-Mch SNela from 
integrated spectral line strengths. However, the timing of 
SNela is important and our current models indicate that 
SNIa yields may contribute to galactic chemical evolution 
earlier than expected. 

Test 6 shows that the effect of removing SNela alto- 
gether is rather large, up to 3.2 x the typical observational 
errors. Clearly this is detectable. The primordial collapse 
models (PCI and PC2) showed smaller differences because 
in those cases SNela do not have much time to make an 
impact before star formation stops. This test illustrates the 
importance of including SNIa yields to correctly interpret 
observed line strengths in galaxies. 

Observe d SNIa rates in galaxies c over about a factor of 
two in range (INavasardvan et alJuOOll . their Table 5). Mod- 
els with double the baseline SNIa contribution were gener- 
ated with proportions given by test 7 in Table [4] Differences 
were found to be up to twice the typical observational uncer- 
tainties so in principle may be detectable for some lines, such 
as iron-sensitive lines. Smaller differences were again found 
for the primordial models. The iron sensitive line-strengths 
are slightly stronger for the doubled SNIa case compared 
with the baseline case (up to 4% of the line-strengths them- 
selves, and not very sensitive to the exact SFH). This is in 
the sense expected, since SNela contribute iron to the cosmic 
cycle. 

Test 8 shows the effect of including a disk wind to 
stabilise accretion. In our models, this is achieved by 
allowing high mass accretion rates to lead to a SNIa, 
in order to test possible SNIa routes first described by 
lHachisu. Kato fc Nomotol (| 19961 ') and subsequent work by 
those authors. As in tests 6 and 7, this shows that the 
line strengths are most strongly affected by the rate of 
SNela. The largest differences were found in models with 
late bursts, as expected from the delayed SNIa contributions 
in these models. 



5 DISCUSSION 

We have shown, using our chemical evolution models, that 
varying parameters that affect the number and type of SNela 
has a greater impact on line strengths than inclusion or ex- 
clusion of yields from other binary stars. This highlights the 
need to focus on the effects of different SNIa models in GCE 
codes. The effects of other binary yields (including common- 
envelope systems that do not become SNIa, CVs, novae etc) 
can largely be ignored, unless specific, low ab undance iso- 
topes , such as 13 C or 15 N are being studied (jlzzard et al.l 
l2006t ). For most GCE applications, where abundances of Fe, 
Mg, C, Ca and overall metallicity have the largest impact 
on line strengths, it is not so necessary to include additional 
binary yields other than from SNela explosions. 

The effects of additional binary yields other than SNela 
explosions are smaller than the observational uncertainties 
on line-strengths. This answers our original question in that 
we do not find that it is necessary to include effects of other 
binary processes besides SNela explosions in GCE models of 
integrated stellar populations in galaxies. This may change if 
a factor of about three improvement in signal-to-noise occurs 
in line-strengths, based on a new, more accurate calibration 
system than the Lick star data. 

We can see from Table [3] that overpopulating with 
SNela (test 8) or excluding SNela altogether (test 6) makes 
the biggest difference to predicted line strengths and these 
are clearly unrealistic and testable. After that, the rate and 
type of SNela are important. The adopted treatment of com- 
mon envelope evolution makes a small difference, but less 
than the observational uncertainties in line strengths. 

iMatteucci et al.l (|2006t ) tested the effects of a more bi- 
modal SNIa rate distribution, using analytical forms, in- 
cluding prompt SNela (< 10 8 yr) and tardy SNela (up to 
10 Gyr) components. They predicted SNIa rates and yields 
for specific SFHs. They found, for ellipticals modelled as 
early, primordial bursts, that the exact SNIa rate distri- 
bution had little effect on SNIa rates at late times. Our 
binary-star yields do include a prompt SNIa (< 10 8 yr) and 
a tardy (continuing) SNIa contribution for both ELDs and 
Mch SNIa types (see Figure |3| . Observational evidence for 
short ( <70 Myr) and long (f ew Gyr) SNIa delay times is 
given bv lAubourg et all ([2008), who look at delay timescales 
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of SNela following star f ormation epochs, in a large sample 
of galaxies. By contrast, iForster et al.l (|2006T ) find that no 
progenitor model can be ruled in or out with current data. 

We find that He WD progenitors for SNela may have 
a small effect in GCE models, however there is as yet no 
observational evidence for these, therefore they are largely 
ignored in this work. If real, their timescales for contribu- 
tions may be longer than other SNela, due to the low mass 
progenitors involved. 

Forcing the baseline SNIa yields to originate from 
Mch COWDs only, by only including scaled-up versions 
of those contributions led to relatively small differences 
in line strengths compared to our baseline model. There- 
fore, whether we model SNela as ELDs or standard Mch 
COWDs makes surprisingly little difference to predicted 
line-strengths in the models now. However, the yields do 
show larger variations with time, which might be detectable 
in resolved stellar populations. Thus for the SNIa explosions 
themselves, GCE modelling of integrated populations does 
not appear to distinguish well between ELDs and standard 
SNela. 

5.1 Effects on average age and composition 
estimates 

Finally, we consider what the implications are for single SSP 
fitting. 

Fitting one SSP to spectral line strengths in order to 
recover luminosity weighted average ages and compositions 
is now a widely used technique in galaxy evolution stud- 
ies. This usually involves a look-up table of SSPs on a grid 
of values in age, metallicity and sometimes abundance ra- 
tio parameter space. The grid spacing we use in our SSP 
fitting routines (based on TMK04 SSPs and code written 
by Robert Proctor) is 0.025 in Log(Age) and Log(Z), with 
coarser grid spacing for abundance ratio. This translates to 
~ 6 % differences between adjacent grid points in age and 
metallicity. Thus the 2% to 7% lower Z in composite popu- 
lations including other binaries (seen in Fig. [3] Z% panels) 
gives rise to one grid spacing lower in luminosity weighted 
average metallicity (Z). 

For luminosity weighted average ages we find no differ- 
ence between estimates from line strengths with or without 
other binaries. The time of the latest starburst (see Fig. |3j 
SFH panels) is approximately recovered in most cases, with 
no systematic differences for fits to line strength with or 
without other binaries. In the case of the model PC2, the 
average age is estimated at 8.4 Gyr in each case, whereas 
the stars in that model really formed between 13.7 and 12.2 
Gyr ago. This highlights the fact that SSP fitting is far more 
strongly affected by aspects other than the effects of other 
binaries besides SNIa. The fact that average age estimates 
are not affected by other binaries is a bit surprising, given 
that all the Balmer indices are stronger when the yields 
from other binaries are included (see Fig. [TJ. However, the 
stronger Balmer indices are offset by other age sensitive fea- 
tures which are systematically weaker (e.g. G4300) and it 
is important to remember that the Balmer indices do have 
some sensitivity to metallicity as well - for a given age, they 
are stronger at lower metallicity. 

Thus average age estimates are unaffected by yields 
from other binaries (such as common envelope systems, CVs 



and Novae) whilst average metallicities are slightly reduced 
(by < 7%) in populations accounting for yields from these 
other binaries. Differences in SNIa types are likely to pro- 
duce larger effects, but SNela are as yet too poorly under- 
stood to address this issue with any accuracy. Future work 
to uncover the progenitors of SNela will help to address this 
large uncertainty in galactic chemical evolution modelling. 
From our present study we can say that other binaries aside 
from SNela have only a small effect on integrated spectra of 
galaxies and so do not normally need to be incorporated to 
interpret galaxy ages and metallicities. 



6 CONCLUSIONS 

In this paper we have investigated the effects of yields from 
binary stars on spectral line strengths in integrated stel- 
lar populations. To do this we combined a single-/binary- 
star population model, which predicts yields as a function 
of time, with a galactic chemical evolution code that self- 
consistently models the spectral line strengths from inte- 
grated populations with different star formation histories. 
Our tested SFHs included stars more than 1.5Gyr old up 
to 13.7Gyr (Table [TJ. Derived line strengths for various as- 
sumed binary-star contributions were compared to a base- 
line model. The resulting differences were compared to typ- 
ical observational uncertainties on line strengths, in order 
to illustrate the relative importance of the various contribu- 
tions and binary effects. We find the following: 

• Populations including other binary star processes as 
well as SNela produce slightly less metal than a population 
consisting of single stars with SNela. This affects magne- 
sium and carbon sensitive spectral features most and iron 
sensitive features least (Figs.[T]and[2]). Hydrogen absorption 
lines are slightly enhanced in populations with other bina- 
ries. However, the differences are all less than 0.9 x typical 
observational uncertainties on line-strengths. 

• Therefore galactic chemical evolution models of unre- 
solved stellar populations, using Lick spectral indices, do not 
need to incorporate yields from binary stars other than SNIa 
explosions because other binary processes (such as CVs, no- 
vae, symbiotic stars and various common envelope binaries) 
do not significantly alter the strengths of observed spectral 
features. 

• Parameters describing common envelope evolution af- 
fect the number of SNela. Varying these parameters (e.g. 
efficiency of orbital energy transfer to the CE) affected ob- 
served line strengths by up to the observational errors in the 
current simulations. Therefore these effects are also quite 
small (less than or equal to observational errors). 

• The exact nature of the explosion (e.g. Sub-Mch or Mch 
COWD) makes less difference than the rate of SNela and the 
time variation of that rate. 

• Reductions in observational errors by about a factor of 
three are needed for line strengths to be sensitive to yields 
from other binary stars besides SNela, or the effects of CE 
evolution, within the bounds of currently acceptable pa- 
rameters. On the other hand, current measurements of line 
strengths are already sensitive to the rates and timescales 
of SNela. However, they are also at least as sensitive to un- 
certainties in massive star evolution. 
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• Luminosity weighted average ages are unaffected by 
yieids from other binaries and iuminosity weighted average 
metaliicities are decreased by iess than 7% due to yieids 
from other binaries. 

This study makes a preliminary investigation into the 
expected contributions of different SNla types to the chem- 
ical evolution of galaxies. The different timescales of SNla 
versus SNII enrichment forms a vital argument in our un- 
derstanding of galaxy evolution from composite stellar pop- 
ulations. Thus, understanding the evolution, timescales and 
chemical contributions of different types of SNela is fun- 
damental for accurate GCE modelling. Therefore it is im- 
portant in future work to further explore the accuracy and 
effects of SNla models and their relative contributions in 
galaxies. Our theoretical yields indicate prompt 10 8 yr) 
as well as longer timescale contributions from SNela (Fig. 2]). 

In this paper we have looked at the effects on spectral 
absorption line-strengths of varying assumptions about the 
contributions of yields from different types of binary stars, 
in chemical evolution models of galaxies. In future work we 
will continue to improve our yield models, guided by ob- 
servational results (e.g. from SuperWasp, CoRoT, Gaia and 
other variability and radial velocity surveys). As more de- 
tailed modelling of individual elements become available we 
will also look at the effects on line strength of individual 
heavy elements such as carbon, which shows the largest dif- 
ferences in line-strengths, due to inclusion of different binary 
star processes other than SNela. 
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APPENDIX A: GCE PARAMETERS 

Parameterisation of the SFH used in our current GCE mod- 
els is aimed at covering as wide a range of types of composite 
galaxy histories as possible with a minimum number of free 
parameters. We also want to specify the latest major burst 
in some detail because more recent star formation has a 
large effect on the luminosity, colours and integrated spec- 
tral properties of the composite stellar population. Since an 
SSP is typically described by three parameters (age, [Fe/H], 
[a/Fe]) we need more than three parameters to describe a 
composite SFH. Empirical fitting of two SSPs typically takes 
six parameters. We aim to improve on that, to describe more 
realistic continuous SFHs. For a GCE model incorporating 
a SFH, we have additional constraints of self-consistency, 
which arise naturally from the model. Therefore the free 



parameters are different from the SSP case. They do not de- 
scribe the attributes of the stellar population directly, but 
instead describe the star formation efficiency and gas flow 
properties as a function of time. 

The early star formation is parameterised by a star for- 
mation efficiency (Co) and a gas flow rate (Fo). To describe 
early star formation, plus details of a late (optional) burst, 
we allow for two interruptions to the history at times T> 0.0. 
These times are represented by Ti and T2, as indicated in 
Fig. IA1I Ti can occur anywhere in the range <Ti < 12.2 
Gyr, i.e. 1.5 Gyr ago or older in the current work, assuming 
an age of the universe of 13.7 Gyr. The duration of the lat- 
est burst is then Di=T2-Ti, as in Table [T] Star formation 
efficiency (Ci) and gas flow rate (Fi) are both set as free 
parameters at Ti. Then at T2 the star formation and gas 
flow are switched off (by fixing C2=0 and F2=0) because 
we are not attempting to model systems with recent or on- 
going star formation. The SSPs used in this current work 
are 1.5 Gyr or older, so this is the limit of how young we 
model stars. This gives us six parameters (the ones given in 
bold above) that can be varied to describe a wide range of 
histories much more realistically than two SSPs (which also 
requires 6 parameters as mentioned above). 

Fig. lAll illustrates the SFH parameters described, along 
with a time sequence covering the age of the Universe. We 
are not sensitive to the exact age of the oldest stars because 
a population that is ~10 Gyr old looks very similar to one 
that is ~14 Gyr old in its spectrum and colours. However we 
are still sensitive to the burst duration in old populations, 
through abundance ratios. Utilizing the six parameters de- 
scribed above we can explore histories ranging from early 
collapse (such as models PCI and PC2 in Table [!}, early 
star formation plus a later burst (such as models LB1, LB2 
and MMW in Table [TJ), through to continuous star forma- 
tion ending recently. The latter models are not explored here 
for modelling effects of binary stars in luminous early-type 
galaxies because [a/Fe] ratio s in these galaxies are n ot con- 
sistent with such models (e.g. iKuntschner et al]|2002l ; PS02) 
and these galaxies do not generally contain significant on- 
going star formation. 



APPENDIX B: LINE STRENGTH 
OBSERVATIONAL UNCERTAINTIES 

We estimate average uncertainties on li ne strengths from ob- 
servations of early- type galaxies (PS02; Dcnicolo et al.ll2005l 
- Dea05 in Table [21 fSanchez-Blazquez et al.ll2006l - S-B06 in 
Table [2} acquired with 4 m class telescopes (e.g. 17 E/S0 
galaxies from PS02, observed with the 4.2 m William Her- 
schel telescope). 

We include the data from PS02 for wavelengths above 
4600A, below this the data are affected by dichroic re- 
sp onse uncertainties. Shor t er tha n 4600A we use the data 
of ISanchez-Blazquez et all (|2006h who provide data on 98 
early-type galaxies as observed with the William Herschel 
telescope or the 3.5 m Calar Alto telescope. 

At wavelengths long er than 5700A errors w ere obtained 
from observations of by iDenicolo etabl (|2005l ). They cata- 
logue 93 early- type galaxies observed with the 2.12 m tele- 
scope in Mexico, so averages of their errors are scaled by 
a factor 0.75 (the ratio of line-strength errors from PS02 
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Figure Al. Schematic illustration of the parameters used to specify a continuous star formation history in our 
galactic chemical evolution model. The six free parameters (without values attached) and three fixed parameters 
(with values) are indicated below the time axis. In Tablc[l]values for the six free parameters are specified to describe 
different models tried in our analysis in this paper. 

to Dea05 in the spectral range 4600 to 5500A). In this way, 
our line-strength errors are estimated from observations with 
the equivalent of a ~ 4m class telescope and a typical ex- 
posure time of 1 to 2 hours. Deeper observations or the use 
of a larger telescope will not drastically reduce these errors, 
because they are dominated by systematic errors in the cali- 
bration to the Lick scale rather t han Poisson errors (e.g. see 
iProctor. Sansom and Reidll2000l. their Table 2; PS02, their 



Table 2; ISanchez-Blazquez et alj |2006| . their Table 4). Ta- 



ble [5] highlights the errors used in the normalisations for 
Fig. El 



